Multi-scale Entanglement Renormalization Ansatz in Two Dimensions: 

Quantum Ising Model 
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We propose a symmetric version of the multi-scale entanglement renormalization Ansatz (MERA) 
in two spatial dimensions (2D) and use this Ansatz to find an unknown ground state of a 2D quantum 
system. Results in the simple 2D quantum Ising model on the 8x8 square lattice are found to be 
very accurate even with the smallest non-trivial truncation parameter. 
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Introduction. — Over the last decade, a rapid de- 
velopment of efRcient methods for simulation of strongly 
correlated quantum systems took place, especially in one 
spatial dimension (ID). It was initiated with the, by now 
classic, paper of White on the density matrix renormal- 
ization group (DMRG) algorithm Recently, the sub- 
ject received new acceleration with the paper of Vidal 01 
who proposed an elegant version of the algorithm based 
on the idea that a state of a ID quantum spin chain can 
be written as a Schmidt decomposition between any two 
parts of the chain. For a generic ground state, but not at 
a quantum critical point, the coefficients of this Schmidt 
decomposition decay exponentially and the decomposi- 
tion can be truncated to a finite number of terms d with 
exponentially small loss of accuracy. The DMRG algo- 
ritms are equivalent to a matrix product state Ansatz [3| 
for the ground state where each spin S is assigned 25* -|- 1 
matrices of size dx d. Each matrix has two indices to be 
contracted with its two nearest neighbors in ID. The ma- 
trix product state can be naturally generalized to two and 
more dimensions by replacing the matrices with higher 
rank tensors to accomodate more nearest neighbors pj. 
These "tensor product states" Q can also be obtained 
as projected entangled pair states (PEPS) [1], the latter 
being a more convenient representation to prove that any 
quantum state can be represented accurately by a PEPS 
for a sufficiently large d. In 2D, unlike in ID, exact cal- 
culation of expectation values in a PEPS is exponentially 
hard with increasing lattice size, but this problem can be 
overcome, at least for open boundary conditions, by an 
efficient approximate method @ which is linear in the 
system size. 

The ability to make efficient and accurate zero tem- 
perature simulations in 2D is of fundamental importance 
for our understanding of strongly correlated 2D quantum 
systems. It is enough to mention possible applications 
to high- Tc superconductors which effectively are 2D sys- 
tems of strongly correlated electrons on a lattice. Their, 
by now classic, Hubbard model [7] has not been solved 
exactly despite staying in the focus of intensive research 
activity for several decades. 

In the context of matrix product states and their gen- 



eralizations, the main difficulty is that all calculations 
are polynomial in the truncation d, but in 2D the degree 
of the polynomial is too high to go far beyond d = 2 
or 3 which, however, may be not accurate enough. A 
possible solution to this problem is the multi-scale en- 
tanglement renormalization Ansatz (MERA) proposed in 
Ref. Q where proper "renormalization" of entanglement 
is shown to reduce the necessary d by orders of magni- 
tude. This economy of the truncation parameter d was 
demonstrated to be truly impressive in the ID quantum 
Ising model where, even at the critical point, a MERA 
with a modest c? = 8 is as accurate as a matrix product 
state with d in the range of a few hundreds 0]. Simu- 
lations with a ID MERA are very efficient because they 
are polynomial in a relatively small d [1, ■ 

The ID MERA is motivated by the following real space 
renormalization group algorithm. Spins on a ID lattice 
can be grouped into a lattice of blocks of two nearest 
neighbor spins. There are two possible choices of block 
chains, A and B, shifted with respect to each other by one 
lattice site. In a decimation step of the renormalization 
group, each A-block is replaced by one effective block spin 
whose Hilbert space is truncated to its d most important 
states. The most important states are the eigenstates 
of an A-block's reduced density matrix with the highest 
eigenvalues. However, to keep d as small as possible but 
without loosing much accuracy before every decimation 
all pairs of nearest neighbor A-blocks are partly disen- 
tangled by 2-spin unitary transformations (disentanglers) 
acting on the 2-spin B-blocks. The disentanglers are op- 
timized to minimize the entropy of entanglement of each 
A-block with the rest of the lattice. They remove entan- 
glement between those pairs of nearest neighbor spins 
which belong to different A-blocks before the A-blocks 
are decimated. The same basic decimation step, includ- 
ing disentanglers, is then applied iteratively to the result- 
ing decimated lattices of block spins. It is worth mention- 
ing, that a similar renormalization group was proposed 
earlier in Ref. [9], but without the disentanglers which 
are essential to keep d reasonably small. 

This renormalization group algorithm motivates the 
MERA. The simplest non-trivial example of a ID MERA 
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is the Ansatz for a periodic lattice of TV = 4 spins: 

"^h^lu Ui^U^, \hk,k,k,) . (1) 

Its graphical representation is shown in panel A of Fig. 
[TJ The repeated indices in Eq. ((T|) imply summation. 
The lowest layer of indices k number basis states of 
the 4 spins. Here the A-blocks are the pairs of spins 
(1,2) and (3,4) and the B-blocks are the pairs (2,3) 
and (4,1). The C/'s are the disentaglers; they are uni- 
tary matrices satisfying unitarity conditions UU'' = 1 

and WU = 1, or C/^^ ('^Sl)* = '5,U3^...4 and 

Kil {uilllY = 4,fc34.fc.. The second layer of in- 
dices j numbers basis states of disentangled spins de- 
fined by e.g. ||j2,j3)) = (^^^^3 1^2, fcs)- The W's are 
isometries or projectors which satisfy orthonormality re- 
lations (W^lfiJ* = Si,u. Their job is to truncate 
the Hilbert space of the disentangled A-block spins to the 
d most important states numbered by the upper indices 
i e {l,...,(i}: for any fixed upper index the matrix 
^hn ^'^^ eigenstate of the A-block reduced den- 
sity matrix in the basis of states of disentangled spins 
||ji, j2))- The eigenstates numbered by indices i become 
basis states |||«i))) = W^jl^^Wji, j2)) of the effective block 
spin. Finally, the top tensor Ti-^^i^, which is normalized 
as Ti^i^ {Ti^i^y = 1, is a quantum state in the basis 
|||u,*2))) of the effective block spins. 

In panel B of Figure [1] we show a generalization of the 
4-spin Ansatz in panel A to a periodic lattice of = 8 
spins. The 8 spins require one more layer of isometries 
and disentanglers. In general, a lattice of A^ = 2" spins 
requires (n — 1) layers of isometries and disentanglers so 
the number of tensors that need to be stored in memory 
is only logarithmic in N . 

In Ref. Q MERA was generalized further to 2D, and in 
Ref. flU it was put in a more general unifying framework. 
In this paper, we propose the alternative 2D Ansatz in 
Fig. [21 In this symmetric Ansatz, 2x2 square plaquettes 
are replaced by effective block spins in each decimation 
step. The symmetric Ansatz is disentangling in a system- 
atic way all those pairs of nearest neighbor (n.n.) spins 
which belong to different 2 x 2-spin decimation blocks, see 
Fig. [3] where the spins on a 2D square lattice are grouped 
into blue and red plaquettes. We propose that in each 
decimation step each blue plaquette is replaced by an ef- 
fective block spin whose Hilbert space is truncated to its 
d most important states, but before each decimation, the 
blue plaquettes are partly disentangled by 4-spin unitary 
disentanglers acting on the red plaquettes. They remove 
entanglement between all those pairs of n.n. spins which 
belong to different blue decimation blocks. Indeed, note 
that in Fig. [3l all links joining such pairs of spins are 
painted red. These red links are naturally grouped into 
red plaquettes and the proposed 4-spin disentanglers re- 
move all the unwanted "red" n.n. entanglement before 
the following decimation. It is essential here that the red 




FIG. 1: In A the ID MERA in Eq. JT]) on a 4-site peri- 
odic lattice of spins. In B the ID MERA is generalized to a 
periodic lattice of 8 spins. 




FIG. 2: The symmetric 2D MERA on a periodic 4x4 lattice. 
The isometries W replace 4-spin square plaquettes with one 
effective block spin in just one decimation step. 



plaquettes are disjoint, because thanks to this all the un- 
wanted "red" entanglement can be removed by the small 
4-spin disentanglers acting on individual red plaquettes. 
Other decimation schemes either do not remove all the 
unwanted n.n. entanglement between different decima- 
tion blocks, or they would require disentanglers acting 
on more than 4 spins. 

The symmetric variant of the renormalization group 
motivates MERA shown in Fig. [2] in the case of 4 x 4 
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FIG. 3: The symmetric decimation in 2D: each blue 4-spin 
square plaquette is replaced by a block spin whose Hilbert 
space is truncated to its d most important states, but be- 
fore this decimation a unitary 4-spin disentangler is applied 
to each red plaquette. The disentanglers remove unwanted 
entanglement between all those (red) nearest neighbor pairs 
of spins which belong to different (blue) decimation blocks. 



periodic lattice. This graph represents the quantum state 



J44J41J14J11 

w*2i . . . 

J24J21J34J31 

rriiiii2 j2i j22 

/Ciifci2fe2lfe22 



Tj331j32]41]42 

k3ik32k4ik42 kssks^k^sk^i' 



IllJl2«2l222 ^ 

J42J43J12J13 

. . X 

J22J23J32J33 

fel3fel4fc23fc24 
rrt33j34j43j44 ^ 



fell 


fcl2 


fcl3 


fci4 


^21 


^22 


^23 


^24 


fcai 


^32 


^33 


^34 


^41 


^42 


^43 


fc44 



(2) 



Here the double subscript indices numerate rows and 
columns of the lattice. A generalization to greater 2" x 2" 
lattices is obtained by adding (n — 2) layers of isometrics 
and disentanglers. 

In this paper we use MERA to find the ground state 



of the spin-i transverse quantum Ising model 



H 



(3) 



on 2 X 2, 4 X 4, and 8x8 periodic square lattices. Here 
X and Z are Pauli matrices. T and all layers of differ- 
ent W and U were optimized to minimize total energy. 
Provided that the minimization preserves all constraints 
on T, U (respectively: normalization, orthonormality 
and unitarity) there is no need to obtain W's as leading 
eigenstates of reduced density matrices and to construct 
U's as disentanglers that minimize the entropy of those 
matrices. The W^s and U's that minimize the energy are 
at the same time good candidates for respectively the 
leading eigenstates and optimal disentanglers. 

In most calculations, we used d = 2 in all tensors i.e. 
the minimal non-trivial value of the truncation parame- 
ter, except for the 8x8 lattice where it was necessary 



to increase the parameter to d = 3, but only in the top 
tensor T near the critical g = 3.04. For any g, the ini- 
tial state for the minimization was the Schrodinger cat 
state I tTT • • • ) + I iii • • ■ ) which is the ground state 
when g ~f 0. This state translates into trivial disenta- 
glers U — 1, the top T having only two non-zero ele- 
ments Tun — T2222 — l/\/2, and all W's being non-zero 
only when VKim — 1^2222 — 1- As we were looking for 
the ground state, we assumed that all tensors T, W, U 
are real. The tensor T and each tensor are quan- 
tum states on a 2 X 2 square plaquette. Each tensor T 
or has four lower indices with each index number- 
ing d states of its corresponding spin. We assume that 
T is symmetric under all exchanges of lower indices that 
correspond to symmetry transformations of the 2x2 pla- 
quette. As each is an eigenstate of a reduced density 
matrix, it must be either symmetric or anti-symmetric 
under each of these symmetry transformations. In all 
considered cases, we found that the lowest energy is ob- 
tained when all W^'s are assumed symmetric under all 
transformations. In this symmetric subspace it is conve- 
nient to parametrize the tensors as (here d — 2) 
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where 

"abed - 



~ means equality up to normalization. Here 
{abcdlv") where the states 
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are a basis of symmetric states on the 2x2 plaquette. 
AaS are imaginary 4-spin hermitian operators invariant 
under the symmetries of the 2x2 plaquette: 

~ Yi+Y2 + Y:, + Yi , 
A2 ~ X1Y2 + Y1X2 + X2Y4 + Y2X4 + 
X3Y4 + Y3X4 + X3Y1 + Y3X1 , 
A21 ~ Y1Y2Y4 + Y2Y4Y3 + Y4Y3Y1 + Y3Y1Y2 . (6) 

Each Aa is a symmetrized sum of tensor products of Pauli 
matrices with each term in the sum including an odd 
number of Y's. They are normalized so that TrAl^A/j — 

Sap- 

The minimized energy is a sum of all bond energies 
Ei,j — {—jgXi — jgXj — ZiZj) . However, thanks to 
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FIG. 4: In A, we compare transversal magnetization (X) on 
the 4x4 lattice obtained from MERA and exact diagonal- 
ization. In B and C, transversal magnetization and nearest 
neighbor ferromagnetic correlator obtained from MERA are 
shown for different lattice sizes. In D, we compare transversal 
magnetization on the 8x8 lattice when d = 2 in all tensors 
and when it is increased to d = 3 in the top tensor with the 
perturbative results from Ref. p^ . 



the assumed symmetry of the tensors, only some of them 
need to be calculated. For example, on the 4x4 lattice 
in Fig. [31 one needs to evaluate only two bond energies: 
^'11, 12 and £'12,13. By symmetry, all other bond energies 
are equal either £11,12 or £12,13 and the total energy is 
= 16£ii,i2 + 16£i2,i3. In a similar way, the 8x8 
square lattice has 6 and, in general, an N x N square 



lattice has ^ + ^ independent bond energies. The to- 



tal number of bonds is 2N'^ so for a large N we save a 
factor of 32 simply by using the assumed tensor symme- 
tries. Thus for a large N, the cost of calculating energy is 
proportional to the lattice size times the cost of calculat- 
ing any bond energy Eij which is logarithmic in N and 
polynomial in d. Here the proof follows similar lines as in 
Ref. [sl . Indices are contracted along causal cones whose 
horizontal cross-section is 3 x 3 (or 4x4) spins when cut 
above (or below) a layer of isometrics W. To avoid the in- 
termediate 4x4 stage we do not apply all isometries first 
and then all disentanglers, but we apply some isometries 
earlier than other gradually including disentanglers i.e. 
we pass through a series of intermediate non-horizontal 
cross-sections never exceeding 11 spins. 

Energy was minimized with respect to the variational 
parameters {iajf^aj?"} in Eq. ^ using different stan- 
dard minimization routines, but the best performance 
was achieved with the simplest steepest descent method 
with gradients of the energy estimated from finite dif- 
ferences. Our calculations demonstrate that energy of 
MERA can be minimized in a fairly straightforward man- 
ner. 

In figure m we summarize our results for 2 x 2, 4 x 4, 



and 8x8 periodic square lattices. Of special interest are 
panels A and D, where we compare transversal magne- 
tization obtained from MERA with exact results on the 
4x4 lattice and perturbative results on the 8x8 lattice. 
On the 4x4 lattice d = 2 was accurate enough, but on 
the 8x8 lattice d in the top tensor had to be increased to 
d = 3. This was necessary because with increasing lattice 
size Ising model develops a critical point at 17 = 3.04 - 
this tendency can be seen in panels B and C. 

In conclusion, we proposed and tested a symmetric 
version of MERA in 2D. Using the smallest non-trivial 
truncation parameter d — 2 in most tensors and fairly 
straightforward optimization methods we obtained sur- 
prisingly accurate numerical results for the ground state 
of the 2D quantum Ising model. This is, we think, an 
encouraging result but, as the Ising model that we con- 
sider is relatively simple, it remains to be seen how well 
MERA can deal with more complicated models. 

Note added. When this paper was in the final stage 
of preparation, the e-print 13] appeared where similar 
symmetric Ansatz was proposed. 
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